Annular magnet system for magnetic resonance spectroscopy

ABSTRACT

A persistent-mode magnet, assembled from superconducting annuli, provides a micro coil NMR, in which compactness and manufacturability are provided for a variety of applications. An annular magnet for micro NMR can include a YBCO-annulus Helmholtz coil, for example, that can energized by a magnet system and then transported for use at a second location with an operating system.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims benefit of U.S. application Ser. No. 11/542,628 filed Oct. 3, 2006 and U.S. Provisional Application No. 60/723,182, filed Oct. 3, 2005, the entire contents of which are incorporated herein by reference.

BACKGROUND OF THE INVENTION

Magnets for magnetic resonance systems and methods for employing such systems in spectroscopy and biomedical applications are in general use. In many such applications it is desirable to provide a system that can produce a magnetic field having the required intensity, distribution and stability for high-resolution analytical spectroscopic applications. Superconducting thin films or wires have been proposed for use in magnets for a variety of applications. However, the cost and field requirements of such systems continues to render achieving the goal of their widespread use in commercial applications still unfulfilled.

There is an ongoing need for improvements in magnet design to provide less expensive magnetic resonance systems having the field strength, spatial distribution and temporal stability needed for a variety of spectroscopic applications.

SUMMARY OF THE INVENTION

The present invention relates to the use of a bulk superconducting material to provide a magnet system for magnetic resonance spectroscopy. In a preferred embodiment, the bulk superconducting material can be in the form of a plurality of annuli that are stacked to provide an annular magnet, the stack being cooled to below the critical temperature of the annuli with a cooling system. The annular assembly surrounds a cavity or bore which can be at ambient or room temperature. The material or sample to be analyzed is placed within the cavity with an NMR probe to obtain magnetic resonance spectra for analysis.

Although there has recently been a great performance improvement in yttrium barium copper oxide (YBCO) bulk disks, for example, such materials have yet to find a useful magnet application. The disk's features, such as improved current densities, provides for a preferred embodiment including a persistent-mode magnet, in which compactness and simple manufacturability are used for high-resolution NMR microspectroscopy or microcoil NMR. In a preferred embodiment an NMR magnet assembled entirely with superconducting annuli is configured to be energized to form a magnetic field at one site, transported, and connected to a magnetic resonance operating system and used at another site. Other superconducting annuli, such as other rare-earth compounds with a current density sufficient for a given application can also be used. The system configuration including critical current densities, the mechanical strength of the superconducting disks, system cryogenics, mechanical integrity, magnet energization procedure, site-to-site transit procedure, stability and protection, temporal field stability, field analysis, and shielding are important elements.

The present invention provides a system operating at 200 MHz or more, preferably at 300 MHz or higher frequencies. The system provides a field of at least 4.70 T (200 MHz) and preferably 7.0 T (300 MHz) or higher. A preferred embodiment, for example, operates at 11.74 T (500 MHz) with an RT bore size between 25 and 60 mm, for example, 43 mm. In this embodiment, liquid neon can be used with an operational temperature of 27.1 K, which can be cooled to a solid 15 K for transport of the annular magnet that has been energized. In a preferred embodiment the annular magnet generates a homogeneous field with a uniformity of at least 0.1 ppm to 0.01 ppm for a 10 mm diameter spherical volume.

A preferred method of using the invention involves the measurement and spectral analysis of a material for use in pharmaceutical and food industries for discovery and development of drugs and foods, and by medical doctors as a in-office outpatient tool for affordable routine clinical diagnostic tests such as blood, other bodily fluid or tissue samples removed from a patient. A material, gas or fluid, can be measured by insertion into the bore of the system with an NMR probe.

DESCRIPTION OF THE DRAWINGS

FIG. 1 is an energizing magnet from an annular superconducting magnet system in accordance with a preferred embodiment of the invention.

FIGS. 2A-2C illustrate an energization procedure in accordance with a preferred embodiment of the invention.

FIGS. 3A-3B illustrate steps of a magnet transport procedure in accordance with the invention.

FIG. 4 is a graphical illustration of H_(z) plots in a preferred embodiment.

FIG. 5 is a B field profile of an annular magnet made in accordance with a preferred embodiment of the invention.

FIG. 6A is a graphical representation of a constant field contour.

FIG. 6B is a graphical representation of a field contour of an annular disk in accordance with the invention.

FIG. 7A is a preferred embodiment of an operating system for the annular magnet system of the invention.

FIG. 7B illustrates a detailed sectional diagram of a magnetic resonance system in accordance with the invention.

FIG. 7C illustrates an embodiment utilizing a coil assembly for magnetic field tuning.

FIGS. 8A-8D illustrate a method of adjusting the magnetic filed intensity in an operating system in accordance with a preferred embodiment of the system.

FIGS. 9A-9D illustrate another sequence for magnetic field tuning.

FIG. 10 illustrates a process sequence of using a preferred embodiment of the invention.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 1 shows a schematic of the system 10 configuration of an annular magnet 16, a cryostat 15 for the annular magnet, an external electromagnet 12 for energizing the magnet; and a coldhead 14. Note that the annular magnet has a room-temperature bore 18 and that the magnet is immersed in a bath of cryogen or surrounded by a volume of solid cryogen. Further details regarding the system can be found in Iwasa et al., “A Persistent Mode Magnet Comprised of YBCO Annuli” Proceedings ASC 2004, Jacksonville, Fla., Oct. 4, 2004, the entire contents of which is incorporated herein by reference. A bulk superconducting material can be a rigid section of material having a cross-sectional area of at least 1 cm², for example, with a thickness of at least 1 mm and a tensile strength as set forth herein. Other rare-earth superconducting compounds of the form Re—Ba—Cu—O such as DyBCO (dysprosium), NdBCO (neodymium), and HoBCO (holmium) which have the desired current density and field uniformity as described herein can also be used.

Critical current density (J_(c)) and tensile strength (σ) data are presented in Table I for a YBCO disk. For a field of 7 T, the material's σ must be at least ˜20 MPa.

TABLE 1 YBCO CURRENT DENSITY & STRENGTH DATA [1] Parameters Value Critical current density, J_(c)A/m² 10⁹ A/m² 77.3 K; 1.08T 0.3 29 K; 17T 3 Tensile strength, σ 293 K [MPa] 53

Most high temperature superconducting (HTS) magnets are operated “dry,” i.e., not immersed in a bath of cryogen but coupled to a cryocooler. However, because the annular magnet in this system can be energized at one site, where a large energizing magnet is stationed, and transported to the user site, in which no such energizing magnet is available, the system, with a fully-energized magnet, must possess a sufficient heat capacity to keep the magnet, during a site-to-site transit, either at a constant temperature in a bath of cryogen or within an acceptable temperature range, e.g., 20-40 K, cooled by a volume of solid nitrogen, for example. The bulk semiconductor material has a current density of at least 10 MA/m² at 77.3 degrees Kelvin and preferably of at least 100 MA/m² at that temperature or below.

Table II presents four heat capacity sources, each with the corresponding magnet operating temperature, T_(op), and a maximum center field, B_(o), achievable. Note that each system, even that relies on a bath of cryogen, is cryocooled, with a coldhead placed at the top of the cryostat. The coldhead either recondenses the cryogen vapor or maintains the temperature of solid nitrogen at 20 K. In the solid nitrogen system, however, T_(op) is permitted to rise up to 40 K during transit and over this temperature excursion the magnet is designed to keep its nominal operating field of B_(o) constant within the system's temporal field stability requirement.

TABLE II HEAT CAPACITY SOURCE; T_(op;) MAXIMUM B_(o) Heat Capacity Source T_(op) [K] Maximum B_(o) [T] 1. Liquid Helium 4.2 >15 2. Liquid Neon 27.1  Up to ~15 3. Liquid Nitrogen 77.3 Up to ~3 4. Solid Nitrogen 20→40* Up to ~1 *Nominal T_(op): 20 K; may rise up to 40 K during transit.

Unlike a conductor-based magnet, in which a reinforcement element may be placed directly in the conductor, co-wound with the conductor, or placed outside the winding, each annulus in this magnet must withstand magnetic stresses. For a field of 10 T, stresses may reach 40 MPa. The YBCO disks for the preferred embodiments of the invention have “trapped” a field of 17 T, or withstood stresses up to ˜120 MPa. For compact microcoil NMR magnets for applications such as drug discovery and development in the pharmaceutical industry, a field of 7.05 T or higher (300 MHz) is deemed adequate; the corresponding stress of 20 MPa are low enough for these annuli. Thus for many applications it is desirable to have a magnetic field generating by the annular magnet at the operating temperature of at least 5 T.

The disks are made using a melt textured growth (MTG) process in which sintered YBCO pellets with the required stoichiometry are partially melted and cooled. The current density at 77 degrees Kelvin or lower is preferably at least 100 MA/m². The material can also be made using liquid phase processing (LPP), quenched and melt growth (QMG) and melt powder melt growth (MPMG) processing. In general, bulk material can be manufacturing using a single crystal growth process from a seed. The resulting product can have a length and width ranging from 20 mm to 100 mm, for example, which is sliced into disks having a thickness of more than 1 mm, typically about 10 mm. Alternatively, in another preferred embodiment, disks can be made using thin annuli of less than 1 mm in thickness that are cut from a thin film or tape of sufficient width and length. Many layers of these thin films are cut into annuli and stacked to provide a single magnet annulus. This can be a laminated structure having in excess of 1000 think film annuli that are each on the order of one micron in thickness.

The procedure for energizing an annular magnet, housed in a cryostat and placed in the bore of an external electromagnet, basically consists of the following steps which are illustrated in FIGS. 2A-2C.

With the YBCO annuli at a temperature above its critical temperature, 93 K, energize the external electromagnet 12 to a level sufficient to enable the annular magnet 16 to trap a field that results in the design center field 20 (FIG. 2A). If shielding steel sheets are to be introduced afterwards for site-to-site transit of an energized magnet, and that can be kept at the user site, the energizing field is preferably adjusted to include the contribution of the steel sheets to the final center field.

While the external magnet is at a constant field, cool the annular magnet to its T_(op) by: immersing it in a bath 22 of liquid helium (T_(op)=4.2 K), liquid neon (27.1 K), or liquid nitrogen (77.3 K); or solidifying liquid nitrogen, already in the cryostat with the cold head, and cooling it further to 40 K, a maximum T_(op) permitted in the system in which the solid nitrogen is the source of heat capacity during the site-to-site transit (FIG. 2B).

Reduce the external field to zero, “field-cooling” each annulus, thereby causing the annular magnet to trap a field 24 (FIG. 2C), which when the magnet is shielded, results in the desired field intensity level with a required spatial distribution. In the solid nitrogen system, further cool the cold body (annular magnet/solid nitrogen) to 20 K, a nominal operating temperature of this system with the cryocooler on. The external magnet is then removed (FIG. 3A) from the annular magnet system for shielding 26 (FIG. 3B) and then packaging.

Though the energizing magnet itself generates a field that is only axially symmetric, satisfying no specific homogeneity requirements, the annular magnet can be configured to create a field with an “NMR-quality” spatial homogeneity.

In the cryogen-cooled system, the volume of cryogen in the cryostat must be sufficient to keep the annular magnet at T_(op) until the time of cryogen refilling at the user site. In the solid nitrogen system, the volume of solid nitrogen must be sufficiently large to keep T_(op) below 40 K. Magnetic field shielding, required for transportation and at the user site, is discussed below.

With HTS magnets, this YBCO annular magnet is stable irrespective of its T_(op), which may range from 4 K to 77 K. For more details regarding YBCO disks for a magnet see Iwasa et al., “Thermal and Magnetic Responses of a Solid Nitrogen/Magnetized YBCO Disk System Undergoing Temperature Cycles in the Range of 8-60 K”, IEEE Trans. Appl. Superconduc. 14 1727(2004) and Tomita et al., “High Temperature Superconductor Bulk Magnets that can Trap Magnetic Fields over 17 Tesla at 29 K”, Nature 421, 517 (2003), the entire contents of these publications being incorporated herein by reference.

An important feature of an annular magnet in terms of protection is that the induced supercurrents in the annuli are not series-connected, rather they, though coupled, can operate essentially independent of one another. Thus, a likelihood of one quenching annulus absorbing the entire stored magnetic energy of the magnet is virtually nill. Quenching may spread out slowly from one annulus to the adjacent annuli, each annulus in the process absorbing as heat at worst slightly more than its own stored energy. Winding materials escape permanent damage if the stored magnetic energy is absorbed as heat by as little as ˜10% of the entire winding.

The induced supercurrent, amplitude and spatial distribution, remains constant as long as each annulus temperature is at or below the field-cooling temperature, T_(fc). To ensure that the condition T_(op)<T_(fc) is maintained, the liquid pressure can be increased when field-cooling: 1.3 atm for helium (T_(fc)=4.5 K); 1.32 atm for neon (T_(fc)=28 K); 1.34 atm for nitrogen (T_(fc)=80 K). Note that for the solid nitrogen cooled system, T_(fc)≦43 K, slightly above the maximum T_(op)=40 allowed during a site-to-site transit.

A field analysis focuses on spatial field homogeneity and a superposition technique. An analytical tool to design a magnet to satisfy spatial field homogeneity is widely available. Along the magnet axis (z), which, in solenoidal magnets, is the direction with the largest change in field, h_(z)≡H_(z)(ç, r=0)/(λ J a₁), where ç≡z/a₁, for a solenoid of winding i.d. 2a₁, o.d. 2a₂, and height 2b and having a constant overall current density; λJ, may be given, with a≡a₂/a₁ and β≡b/a₁, h _(z)(ç)=F(a,β)[1+E ₂(a,β)ç² + . . . +E _(2j)(a,β)ç^(2j)]  (1) For a solenoid its center field, H_(z)(0, 0), may be given, with F(α,β)≡1n[(α+√{square root over (α²+β²)})/(1+√{square root over (1+β²)})], by H _(z)(0, 0)=λ J a ₁ F(a,β)   (2) Coefficients, E₂(a,β), . . . , E_(2j)(a,β), are known, e.g.,

$\begin{matrix} {{E_{2}\left( {\alpha,\beta} \right)} = {\frac{1}{2\beta\;{F\left( {\alpha,\beta} \right)}}\left\lbrack {\frac{1}{\left( {1 + \beta^{2}} \right)^{1.5}} - \frac{\alpha^{3}}{\left( {\alpha^{2} + \beta^{2}} \right)^{1.5}}} \right\rbrack}} & (3) \end{matrix}$ A superposition technique may be applied to perform a field analysis on the n^(th) pair of annuli symmetric about z=O. To account for a spatial variation of supercurrent density, J_(s)(z, r, B), within each of the n^(th) annulus pair, it is possible to segment each into m sub-annuli, within which J_(s)(z,r,B) is uniform, i.e., J_(s) ^(m). A pair of k^(th) sub-annuli may in turn be represented as a sum of two uniform-current density sub-solenoids k and k′ of opposite current directions. Note that the size of each uniform-J_(s) ^(m) sub-annulus depends on the extent of spatial variation of J_(s)(z, r, B) within each annulus, which in turn must be computed by numerical analysis.

Two dimensionless fields, h_(z)(ç)|_(k) ^(n)≡H_(z)(ç)|_(k) ^(n)/(λJ_(s) ^(m)a₁)_(k) ^(n) and h′_(z)(ç)|_(k) ^(n)≡H′_(z)(ç)|_(k) ^(n)/(λJ_(s) ^(m)a_(a))_(k) ^(n), of sub-solenoids k_(n) and k_(n)′ in the nth annulus pair, respectively, are given by:

$\begin{matrix} {{{h_{z}(ϛ)}❘_{k}^{n}} = {{F\left( {a_{k}^{n},\beta_{k}^{n}} \right)}\left\lbrack {1 + {\sum\limits_{j = 1}^{\infty}{{E_{2j}\left( {a_{k}^{n},\beta_{k}^{n}} \right)}ϛ^{2j}}}} \right\rbrack}} & \left( {4a} \right) \\ {{{h_{z}^{\prime}(ϛ)}❘_{k^{\prime}}^{n}} = {{F\left( {a_{k}^{n},\beta_{k^{\prime}}^{n}} \right)}\left\lbrack {1 + {\sum\limits_{j = 1}^{\infty}{{E_{2j}\left( {a_{k}^{n},\beta_{k^{\prime}}^{n}} \right)}ϛ^{2j}}}} \right\rbrack}} & \left( {4b} \right) \end{matrix}$

where a_(k) ^(n)=(a₂/a₁); β_(k) ^(n)=(b/a₁); β_(k′) ^(n)=(b′/a₁). H_(z)(0,0)|_(k) ^(n) and H′_(z)(0,0)|_(k′) ^(n) for respective sub-solenoids are given by: H _(z)(0,0)|_(k) ^(n)=(λJ _(s) ^(m) a ₁)_(k) ^(n) F(a _(k) ^(n),β_(k) ^(n))   (5a) H′ _(z)(0,0|_(k′) ^(n)=(λJ _(s) ^(m) a ₁)_(k) ^(n) F′(a _(k) ^(n),β_(k′) ^(n))   (5b)

The combined dimensionless z-axis field of pair of k^(th) sub-annuli in a pair of n^(th) annuli, [h_(z)(ç)]|_(k) ^(n), is given by: [h _(z)(ç)]|_(k) ^(n) =h _(z)(ç)|_(k) ^(n) −h′ _(z)(ç)|_(k′) ^(n)   (6) From (4) and (6) we have an expression for [h_(zk)(ç)]_(k) ^(n):

$\begin{matrix} {{\left\lbrack {h_{zk}(ϛ)} \right\rbrack ❘_{k}^{n}} = {\left\lbrack {{F\left( {a_{k}^{n},\beta_{k}^{n}} \right)} - {F^{\prime}\left( {a_{k}^{n},\beta_{k^{\prime}}^{n}} \right)}} \right\rbrack + {\sum\limits_{j = 1}^{\infty}{\left\lbrack {{{F\left( {a_{k}^{n},\beta_{k}^{n}} \right)}{E_{2j}\left( {a_{k}^{n},\beta_{k}^{n}} \right)}} - {{F^{\prime}\left( {a_{k}^{n},\beta_{k^{\prime}}^{n}} \right)}{E_{2j}^{\prime}\left( {a_{k}^{n},\beta_{k^{\prime}}^{n}} \right)}}} \right\rbrack ϛ^{2j}}}}} & (7) \end{matrix}$

The z-axis field of a magnet comprising n^(th)-annulus pairs can be computed by using equation (7) and summing up each field contribution: first uniform-J_(s) m sub-annulus pairs within the n^(th)-annulus pairs, then n annulus pairs. Adjustable design parameters to meet a set of spatial field homogeneity requirements include inner diameter (i.d.), outer diameter (o.d.), and thickness of each annulus (in pair); spacing between adjacent annuli; and number of annuli in the magnet. As discussed below, each annulus-pair must be substantially identical not only in size but also, more importantly, J_(s)(z, r, B), a set of demanding requirements but not unlike that of strict manufacturing specifications of high-resolution NMR magnets.

The fringing field of an energized annular magnet outside the package housing it for shipment is preferably “small”. For a viable compact system, its fringing field must be ≦5 gauss for the user at ˜1 m from the magnet center.

In this embodiment, because of a coldhead at the top end of the cryostat, shielding steel sheets cannot be placed at the top or bottom end of the system—the system requires a symmetrical placement of magnetic materials. Thus steel annular sheets may be placed outside the cryostat either: before energization (option 1); or after energization but before shipment (option 2) to the user site. Option 1 requires an external magnet with a large room-temperature bore; because the system is already energized. Option 2 requires a careful placement procedures of the steel sheets.

In a first embodiment, the following describes a field analysis and measurements on a two-annulus Helmholtz coil. Each annulus, 25-mm i.d.; 48-mm o.d.; and 5-mm thick, was prepared from a 48-mm o.d. YBCO rod. Although, as stated above, each annulus pair is preferably identical, the two annuli used in this measurement had J_(s)(z,r,B) distributions that differed slightly from each other. Such field disparities can be corrected using a magnetic field tuning device and method described herein.

In an example of the system is shown in FIG. 1, a two-annulus Helmholtz coil, with an adjustable gap distance, is used in place of an annular magnet. Here, an energizing magnet, 140-mm id., 320-mm o.d., 188-mm height, is wound with copper conductor and operated in a bath of liquid nitrogen. The Helmholtz coil's H_(z)(z) is measured by a motor-driven search coil, at the central room-temperature bore and moved in the z-direction at a constant speed.

First, immerse the copper magnet in liquid nitrogen and, while the YBCO-annuli are still dry and above 93 K, energize the copper magnet to a constant field. Next, immerse the Helmholtz coil in liquid nitrogen, turning the YBCO annuli superconducting. Then, discharge the copper magnet, field-cooling the YBCO annuli. Then, move the search coil along the z-axis with a constant speed and integrate, digitally, the search coil output signal to generate an H_(z)(z) plot, which is recorded.

In these measurements, comparison between experiment and analysis was limited to H_(z)(z). Thus, the analysis did not invoke equation (7); rather each annulus was segmented into sub-annuli to account for a spatial variation in the induced supercurrent distribution, J_(s)(x,r,B), and the field contribution from each sub-annulus was summed up along the z-axis to compute H_(z)(z).

FIG. 4 shows a set of three H_(z)(z)plots, measured (solid), computed (dashed), and “weighted” (chain dotted), field-cooled at 0.08 T and 77.3 K. Here the measured plot corresponds to H_(z)(z) averaged over the search coil active volume; the computed plot represents H_(z)(z) for each axial point; and the weighted represents H_(z)(z) averaged over the search coil active volume, axial as well as radial.

Although agreement between measured and weighted plots is excellent, the measured H_(z)(z) is not completely symmetric about z=0. This field asymmetry is caused by a difference, however small, of J_(s)(z,r,B) distributions of the two annuli. Thus, to synthesize an NMR-quality magnetic field with an annular magnet, induced by an external magnetic field of a nominal axial symmetry, each annulus-pair, though not all annuli, preferably have identical J_(s)(z,r,B) (i.e., less than 1% variation in current density). For this Helmholtz coil this field asymmetry, as may be expected, became more pronounced as the energizing field was increased and the J_(s)(z,r,B) difference became greater.

These measurements indicate that for magnets where compactness and simple manufacturability are vital, such as for microcoil NMR, this type of “persistent-mode” magnet, assembled entirely from YBCO annuli is suitable for this application. This analysis has shown that to synthesize a field with an NMR-grade spatial homogeneity, each annulus pair is preferably substantially identical not only in size but more crucially J_(s)(z,r,B).

FIG. 5. Shows a contour plot of B of a field cooled system at 77.3 K, for example. FIG. 6A shows constant field contours resulting from the data of FIG. 5. The contours in FIG. 6A indicate that the supercurrent distribution at each radial location in this annulus is not angle independent, which is preferably true in 300 MHz to 600 MHz systems. FIG. 6B illustrates the field contours of an annular disk in accordance with the invention having further improvements in field intensity and uniformity.

In an example, consider that the average critical current density {tilde over (J)}_(c)(77.3 K) is related to the annulus center flux density, B_(z)(0,0,0), which is derivable from B_(z)(0,0,b). Thus:

$\begin{matrix} {{B_{z}\left( {0,0,0} \right)} = {\mu_{o}{{\overset{\sim}{J}}_{c}\left( {77.3\mspace{14mu} K} \right)}a_{1}\beta\;{\ln\left\lbrack \frac{\alpha + \sqrt{\alpha^{2} + \beta^{2}}}{1 + \sqrt{1 + \beta^{2}}} \right\rbrack}}} & (8) \end{matrix}$ For a uniform current density solenoid (or annulus), B_(z)(0,0,b) and B_(z)(0,0,0) are related: B _(z)(0,0,b)=½B _(z)(0,0,0;β′)   (9) where B_(z)(0,0,0;β′) is the center field of a solenoid (or annulus) of the same {tilde over (J)}_(c)(77.3 K) except with β′=2β. Equation 9 may be derived from simple superposition.

Combining Eqs. 8 and 9, and substituting 2β for β in Eq. 8, we obtain:

$\begin{matrix} {{B_{z}\left( {0,0,b} \right)} = {\frac{1}{2}\mu_{o}{{\overset{\sim}{J}}_{c}\left( {77.3\mspace{14mu} K} \right)}{a_{1}\left( {2\beta} \right)}{\ln\left\lbrack \frac{\alpha + \sqrt{\alpha^{2} + \left( {2\beta} \right)^{2}}}{1 + \sqrt{1 + \left( {2\beta} \right)^{2}}} \right\rbrack}}} & (10) \end{matrix}$ Solving Eq. 10, for {tilde over (J)}_(c)(77.3 K), with α=1.84, β=0.2, and a₁=0.0125 m, the critical current density is:

$\begin{matrix} {{{\overset{\sim}{J}}_{c}\left( {77.3\mspace{14mu} K} \right)} = {\frac{B_{z}\left( {0,0,b} \right)}{\left( {4\pi \times 10^{- 7}\mspace{14mu}{H/m}} \right)\left( {0.0125m} \right)(0.2)} \times}} \\ {\frac{1}{{In}\left\lbrack \frac{1.84 + \sqrt{(1.84)^{2} + {4(0.2)^{2}}}}{1 + \sqrt{1 + {4(0.2)^{2}}}} \right\rbrack}} \\ {= {545.4{B_{z}\left( {0,0,b} \right)}\mspace{14mu}{{MA}/m^{2}}}} \end{matrix}$

A system 40 for operation at 500 MHz is shown in FIG. 7A. In this embodiment there are 24-28 annuli 42, for example, which are each separated by a spacer 44 such as a high strength thermally conductive metal such as a copper-silver alloy. The system has steel annuli sheets 46 placed outside the cryostat for field shielding. The shield can be thermally anchored to the first stage of the coldhead. The system uses metal coils 48 such as room temperature copper winding to adjust the field prior to use. A more detailed illustration of a system 60 in accordance with a preferred embodiment of the invention is shown in FIG. 7B.

FIG. 7B shows a schematic drawing of the cross-sectional view of the preferred embodiment of the present invention at the user site. This desk-top NMR magnet system 60 can consist of 16 annuli 62, each 55-mm i.d. and 10-mm thick. A number of annulus pairs, at the upper and lower sections, have an o.d. of 80 mm, while the remaining pairs of annuli in the center section have o.d.'s less than 80 mm, creating a “notch” effect. As schematically indicated in FIG. 7B, the entire magnet assembly is clamped axially with a support assembly 64 such as tie rods placed outside the annuli. The system can include thermal fins 80 and radiation shield 86.

Where the system has a uniform current density, 2a₁=55 mm, 2a₂=80 mm, and 2b≅187 mm (here the spacers are included), i.e., a≅1.455 and β≅3.364, computing J requires the annuli for the system to generate a central field, B_(z)(0,0,0), of 4.7 T: λ J≅320 MA/m². Because in an NMR-field-quality annulus magnet, B_(z)(0,0,0)=4.7 T is nearly equal to the “peak” field, 320 MA/m² is the average critical current density required in each annulus of the prototype, field-cooled at 27.1 K and operated nominally at 15 K. Note that with 320 MA/m² is ˜ 1/10 of {tilde over (J)}_(c)≅3000 MA/m² at 29 K, can be achieved with the YBCO disks.

A spacer is placed between adjacent annuli to satisfy three requirements: 1) field homogeneity; 2) mechanical reinforcement; and 3) thermal stability. A high-strength, conductive sheet of Cu/Ag alloy can be used for the spacer material.

At 4.7 T the maximum Lorentz-force induced tensile stress in each annulus, 80-mm o.d., 55-mm i.d., and 10-mm thick, is 45 MPa, well within the capacity of YBCO disks. High-strength (ultimate strength above 500 MPa and preferably close to 1,000 MPa; Young's modulus ˜150 GPa) spacers takes some radial load off the annuli. The magnet assembly can be precompressed with tie rods to supplement the axial magnetic clamping in the magnet, which enhances the friction between the annuli and the spacers, helping in turn to ensure that the spacers share the load.

During the field-cool process when the energizing electromagnet is discharged, AC losses, primarily hysteretic, are generated in the annuli. Fortunately, the field discharge rate can be as slow as required to limit the dissipation-induced temperature rise, despite the low thermal conductivity of YBCO, in the annuli. The hysteretic energy density e_(hy) [J/m³], may be given by:

$\begin{matrix} {e_{hy} = {\frac{1}{V_{o}}{\int_{V}{\left\lbrack {\int_{0}^{\tau_{d}}{{\overset{->}{E} \cdot \overset{->}{J}}{\mathbb{d}t}}} \right\rbrack{\mathbb{d}V}}}}} & (11) \end{matrix}$ where {right arrow over (E)} is the electric field induced by the time-varying magnetic flux; {right arrow over (J)}, here |{right arrow over (J)}|=J_(c), is the current density in the annulus; V_(o)[=2π(a₂ ²−a₁ ²)b] is the annulus volume; and τ_(d), is the discharge time. Assume that J_(c), (B) is constant at 1000 MA/m². (Note that an average current density of only 320 MA/m² over the entire cross section of each annulus is required to produce 4.7 T; use of 1000 MA/m² is appropriate for current confined only to ˜⅓ of the annulus volume.) With (a₂−a₁)=12.5×10⁻³ m and B_(fc)=5 T, obtain: e_(hy) 40×10⁶ J/m³ or 40 J/cm³: the average heat generation density is q_(hy)=e_(hy)/τ_(d) For an annulus of thickness 2b (10 mm), assume that heat flows only axially (upward and downward) towards the nearest Cu—Ag spacers. Then the maximum temperature rise, ΔT mx, occurs at the midplane of each annulus, given by:

$\begin{matrix} {{\Delta\; T_{mx}} = {{\frac{q_{hy}}{2k}b^{2}} = {\frac{e_{hy}}{2k_{an}\tau_{d}}b^{2}}}} & (12) \end{matrix}$ where k_(an) is the thermal conductivity of YBCO. With ΔTmx=0.5 K, k_(an)=1.5 mW/cm K (“epoxy” at 20 K), b=0.5 cm, and e_(hy)=40 J/cm³, we obtain: τ_(d)=6600 s; a discharge time of ˜2 hours should be long enough to keep each annulus within 0.5 K of 27.1 K. Each Cu—Ag (k˜10 W/cm K) spacer of ˜1-mm thickness can easily conduct this heat to the liquid neon that surrounds the annulus magnet during the field-cooling process.

Because virtually the entire axial flux of the bore flows through the steel assembly of D_(so) o.d. and D_(si) i.d., conservation of flux results in:

$\begin{matrix} {{\frac{\pi}{3}\left( {a_{2}^{2} + {a_{2}a_{1}} + a_{1}^{2}} \right)B_{o}} \simeq {\frac{\pi}{4}\left( {D_{so}^{2} - D_{si}^{2}} \right)\mu_{o}M_{s}}} & (13) \end{matrix}$ where B_(o) is the central field. The left-hand side of Eq. 13 represents the total flux in the magnet bore (π B_(o) a₁ ²) plus that in the annular space between a₁ and a₂, where B(r)=B_(o)(a₂−r)/(a₂−a₁), decreasing linearly from B_(o) a₁ to 0 at a₂; M_(s) is the magnetization of steel, low enough to make its effective μ≧100 μ_(o). For a preferred embodiment, (B_(o)=4.7 T; 2a₁=55 mm; 2a₂=80 mm), with μ_(o) M_(s)=1 T: D_(si)=229 mm and D_(so)=272 mm. A total of 150 these annuli, each 2-mm thick, for a mass of 40 kg, are used.

Steel annuli can be placed one at a time, manually. The steel assembly 88 can be to 300-mm long, nearly twice as long as the annulus magnet 62: unsupported, the midplane of the steel assembly is about 40 mm lower than that of the annulus magnet. The tie-rod assembly 64 attached to the bottom of the cryostat's upper section, supports the steel assembly. Conical spring washers maintain the clamping force constant despite difference in thermal contraction between system components.

In the r-direction, the coupling force between the two assemblies is unstable. The destabilizing radial force is zero when the two are symmetric in the φ and z-directions, but the force increases with radial displacement. If both bodies are symmetric in φ, a 1 mm offset sideways generates a lateral force of about 200 N. Because field homogeneity itself requires concentricity, it can be met by tight dimensioning and structural rigidity of the annulus magnet, the cryostat, and the steel assembly. This concentricity requirement is no different from that for a standard NMR magnet system, in which one such example is between the main magnet and its field shielding elements, a set of active coils or ferromagnetic tiles.

The coldhead 66 has two stages: the first stage 78 (30 K) for the radiation shields and the second stage 75 (15 K) for the magnet and Ne 84. In this embodiment, the coldhead remains with the cryostat even when it is shut off. To minimize heat input from the shutoff coldhead, each coldhead stage can be disengaged from its cooling target by a thermal switch 65 with actuator 68. Also, to minimize its vibration on measurement, the coldhead is mounted on the cryostat top plate through a vibration isolator 76. Of course, to provide a noise-free measurement environment, at least for up to two days, the coldhead may be shut off. A flexible thermal path 82 couples stage 75 to switch 65. A charcoal getter 77 is also used. Insulation 74, radiation shield 72 and a fill/vent port 70 for the coolant are also used.

The cryostat has a room-temperature bore of 38 mm. As shown in FIG. 7B, it has a space large enough to house the coldhead in its upper section. The tail section uses standard cryostat construction techniques to ensure concentricity of the tubes. The cryostat will house ˜4 liters of solid Ne, so that it can take 48 hours for the entire solid Ne volume to warm up from 15 K to 24 K with a net average heat input of 0.3 W. The system can stay within this heat input limit, provided the radiation shields, disengaged from the 1^(st) stage of the shut-off coldhead, remain below 60 K.

When the coldhead is disengaged, heat leaks raise the temperature of the shields from 30 K. The total radiation heat input to the shields of a total surface area of 6,600 cm² facing room temperature is 2 W over the temperature range 30-60 K. The conductive heat input to the shields through the support tubes (thin-walled stainless steel) over the same temperature is 1 W. The convective input is negligible, because the cryostat vacuum is kept below 10⁻⁵ torr with a charcoal getter placed in the cryostat vacuum space.

A heat input of 3 W increases the 20-kg copper shields' thermal energy by 520 kJ over a 48-h period, of which 50 kJ is absorbed by the shields. The large enthalpy density, 50 J/cm³, of solid nitrogen from 30 K to 60 K (63.1 K melting point), absorbs the rest (470 kJ): a 10-liter volume solid nitrogen (0.9 kg) is attached to the radiation shields.

The heat input (radiation and conductive) to the cold body consisting of the annulus magnet and a volume of solid neon is on average 0.3 W as the radiation shields warm up from 30 K to 60 K over a 48-h period. The required volume of solid neon to absorb a total enthalpy input of 52 kJ as it warms up from 15 K to 24 K is ˜4 liters.

When neon freezes at 24.57 K, it contracts 15.9% by volume, or ˜5% linear. From 24.57K down to 15 K, it contracts another 2.1% by volume, or ˜0.7% linear. Similarly, nitrogen upon freezing at 63.1 K contracts 2% by volume, or ˜0.7% linear, and between 63 K and 15 K, contracts another 17% by volume, or ˜5% linear. Since there have been no adverse contraction (or expansion) effects, structurally or thermally, of solid N₂ systems built and operated in the range 10-77.3 K, there are no undesirable effects, structurally or thermally, in any solid Ne-cooled magnet systems. An RF coil probe 90 and sample can be inserted into the bore for spectral measurements. Gradient field coils can also be included in the probe or integrated into the magnet system. A spinner can also be included in the probe to rotate the sample container which can contain a gas or liquid. The coils and magnet system can be controlled by an operating system 100 such as a computer. The system 100 can include a data processor for data storage, processing of spectral data and display. The system 100 can be programmed to perform a pulse sequence excitation of the sample.

FIG. 7C shows a detailed view of an annular magnet system 120 having a coil 128 with windings surrounding each annulus 140. The coil 128 provides for “tweaking” or magnetic field tuning of axial field variations. The coil is surrounded by shielding plates 142. The coolant 126 such as solid neon is housed in cryostat 124 which surrounds the room temperature bore 122. Shim coils can also be used with the system to make adjustments to the magnetic field distribution and can be placed in the annular cavity or bore about the sample to be measured. These can include x,y and xy corrections. Ferromagnetic tiles can also be placed in the annular cavity to reduce x,y and error harmonics.

In the sequence illustrated in FIGS. 8A-8D a process for adjusting the field H(r) is illustrated. The shaded area is the cross-section of the annulus and the center is the annulus bore. In FIG. 8A, the dashed line shows the nonsuperconducting annuli in the external field of the energizing magnet and the solid line after field cooling. FIG. 8B shows the field after application of a positive field and a return to zero thereof. This indicates that a small positive adjustment in the field has no effect on the field in the bore. FIG. 8C shows a dashed line illustrating the application of a negative adjusting field and FIG. 8D shows the field after the adjusting field is returned to zero indicating a permanent adjustment in the field within the bore. Each annulus can have separate coils to provide a substantially separate adjustment for each annulus.

In another embodiment of this field-tuning process is a correct-and-set type method. FIG. 9A-9D shows four sets of H_(z)(r) plots within an annulus of a “long” annulus assembly. Each plot is based on the critical-state model for a “Bean slab” of width 2a.

The shaded areas represent the cross section of an annulus (2a₁; 2a₂), and the white center region the annulus bore of diameter 2a₁—as stated above, it is an infinitely long and deep slab of width 2a₂, with the central part of width 2a₁ removed: each annulus is in the critical-state within a radial penetration depth of δ_(p), indicated in FIG. 9A, where the current density, J_(c), is induced. In FIG. 9A the dotted lines correspond to the field profile when the entire annulus is in the critical state. As noted above, for the embodiment where the annulus magnet generates a central field of 4.7 T, it requires an average current density of 320 MA/m² uniformly distributed over the entire cross section of each annulus.

As in each H_(z)(r) profile shown here, dH_(z)(r)/dr=J_(c), in which J_(c) is assumed field-independent here. In the actual case, because J_(c) is field-dependent, the profile, and hence the bore field distribution, must be computed numerically.

The dashed horizontal line represents H_(z)(r) of the normal-state annulus in the field of the energizing electromagnet. The solid lines represent H_(z)(r) of the field-cooled annulus. (In reality there is virtually no difference in field amplitude between that represented by the dashed line and that by the solid horizontal line.) Realistically, the annulus is in the critical state only within the penetration depth of δ_(p). Because J_(c)=10⁹ A/m², and δ_(p)=H_(o)/_(Jc) (where H_(o)=B_(o)/μ_(o)), δ_(p)≅3.7 mm<<12.5 mm=(a₂−a₁), for B_(o)=4.7 T—for clarity, δ_(p) is drawn greater than its actual value.

The profile in FIG. 9B corresponds to the annulus after it is subjected to a small positive tweaking field, H_(tf). Its amplitude is drawn much greater than its actual size. Within a layer of δ_(p|)=H_(tf)./J_(c), a current is induced in the direction opposite to that already induced in the annulus, shielding the interior. Its contribution to the central field will be slightly reduced, because the net current flowing in the annulus that generates the central field is reduced by the presence of the current within this layer that flows in the opposite direction. When the tweaking field is reduced to zero, H_(z)(r) returns to the original profile (the solid lines in FIG. 9A).

When a small negative tweaking field is applied as seen in FIG. 9C, the field, as indicated by the dashed lines, penetrates further into the annulus. Its contribution to the central field will be slightly enhanced because the net induced current in the annulus is increased.

When the negative tweaking field is returned to zero, as in FIG. 9D the current is induced at the outermost layer of δ_(p−) in in the direction opposite to that of the original current. The original current layer (δ_(p)) is thus shifted inward by and and a new current layer δ_(p−) is created: a net decrease in the central field.

The field-tweaking coil for the prototype (16 annuli) magnet consists of 16 copper-wound coils, each of which can be energized independently. The tweaking procedure can be iterative and make “gross” correction to the first few axial terms, i.e., z, z², and z³. Further refinement must be dealt with by the combination of room-temperature shim coils and ferromagnetic shims placed in the annular space in the magnet bore. With δ_(p−) given by δ_(p−)=H_(tf)/J_(c), we compute δ_(p−) to be 0.08 μm for μ_(o) H_(tf)=10⁻⁴ T and J_(c)=10⁹ A/m².

Before the application of a tweaking field, the current element in the annulus has an i.d. of 2(a₂−δ_(p)), an o.d. of a₂, and a length 2b (not to be confused with an annulus thickness of 2b). When the current layer shifts inward by δ_(p−)<<a₁, the corresponding dimensions become, respectively, 2(a₂−δ_(p)−δ_(p−)), 2(a₂−δ_(p−)), and 2b (same length). When a current element's a₁, a₂ and 2b are changed to a′₁=a₁−δ and a′₂=a₂−δ, the new current element parameters become: α′=(α−ε)/(1−ε) and β′=β/(1−ε), where ε=δ/a₁. For ε<<1, we may derive ΔB_(o), the change in the central field, from Eq. C1 with B_(z)(0,0,0)=B_(o):

$\begin{matrix} {{\Delta\; B_{o}} = {\left( \frac{\delta_{p^{-}}}{a_{1}} \right){B_{o}\left\lbrack {\frac{1}{\sqrt{1 + \beta^{2}}} - \frac{1}{\sqrt{\alpha^{2} + \beta^{2}}}} \right\rbrack}}} & (14) \end{matrix}$ Note that, as expected from Eq. 8, when a₁ and a₂ are decreased, the central field increases. For the original current element dimensions of a₁=40 mm−3.7 mm=36.3 mm (J_(c)=10⁹ A/m2), a₂=40 mm, and 2b=185 mm, with δ=δ_(p−)=0.08 μm, B_(o)=4.7 T, α=1.102 and β=2.548, one obtains: Δ B_(o)≅5×10⁻⁴ gauss, a negligible change.

To this increase in the central field, which, here at 5×10⁻⁴ gauss, is really negligible, there is a subtracting flux density, B_(o−), generated by the newly created thin current layer of δ_(p−)<<a₁ flowing in the opposite direction at the rim, given approximately by: B_(o)−≃μ_(o)J_(c)δ_(p)−  (15) which for this case is equal to 1 gauss. That is, the net decrease in the central field after application of B_(ft)=1 gauss in the negative direction is 1 gauss.

One field-tuning coil, a single turn of a 1-mmø copper wire wrapped around each annulus, generating 1 gauss, requires 3.2 A. At 15 K(p_(cu)=2×10⁻⁸ Ωcm), Joule dissipation for this 1-gauss copper coil is 0.7 mW, a negligible level. (Nevertheless, the cryocooler will be kept running during field-tweaking). Because 16 pairs of 1-mmø copper conduct a lot of heat to the cold anchor (1st stage radiation shield), they will be disconnected from the anchor when the coldhead is disengaged for shipment or quiet measurement environment.

A field of no greater than 10⁻⁴ T (1 gauss) can be used to make the axial field NMR-quality. To evaluate this tuning-field coil, apply a tuning field to the annulus magnet, i.e., superconducting at 15 K but not yet field-cooled. In this state, the field distributions would be different from those shown in FIGS. 9A-9D, but each on-and-off tuning-field sequence of ˜1 gauss will result in a trapped field of the same order of magnitude, which, in the absence of a background field of 4.7 T, are detectable with a Hall probe.

This measurement with the field-tuning coil can also quantify proximity effects on the neighboring annuli by the tuning field of the targeted annulus. Because the winding diameter of each 1-gauss field-tuning coil is 81 mm (except for those annuli in the center zone), it generates axial fields of 0.90, 0.68, 0.47, and 0.31 gauss, respectively, at the centers of successively more distant annuli. Clearly, each tuning field targeted at a particular annulus will likely affect quite significantly the magnetization of at least its two nearest annuli.

A method of using the magnetic resonance system in accordance with the invention is illustrated in the process sequence 200 of FIG. 10. After the annular magnet is energized 202 the magnet is optionally transported and connected 204 to the operating system of a magnetic resonance system. A sample is obtained such as a pharmacological agent or drug, or a bodily sample from a mammalian animal or patient. The sample is inserted 206 into the system and a spectrum is obtained 208. The system can then analyze the spectrum to determine a characteristic 210 such as a diagnosis of a condition of the patient or the suitability of a drug to treat a particular condition.

While the invention has been described in connection with specific methods and apparatus, those skilled in connection in the art will recognize other equivalents to the specific embodiments herein. It is to be understood that the description is by way of example and not as a limitation to the scope of the invention and these equivalents are intended to be encompassed by the claims set forth below. 

1. A magnetic resonance system comprising: a a persistent mode magnet system having a plurality of superconducting annular layers, the annular layers comprising a rare-earth barium copper oxide compound having a tensile strength of at least 20 MPa and a current density of at least 10 MA/m² in a temperature range of 4k to 80k; a cooling system that cools the plurality of superconducting annular layers to an operating temperature such that the magnet system generates an operating magnetic field of at least 4.7T in a system cavity; and an energization system that energizes the persistent mode magnet system.
 2. The system of claim 1 further comprising an operating system having a conductive magnetic field tuning coil.
 3. The system of claim 2 wherein the conductive coil comprises a winding around each of the plurality of layers.
 4. The system of claim 2 wherein the conductive coil adjusts a magnetic field generated by the persistent mode magnet system.
 5. The system of claim 1 wherein system operates at a frequency of at least 300 MHz.
 6. The system of claim 1 wherein the cavity has a diameter of 25 to 60 mm.
 7. The system of claim 1 wherein the cooling system comprises a cryostat and a coolant.
 8. The system of claim 7 wherein the coolant comprises a liquid or solid cryogen.
 9. The system of claim 1 wherein the cooling system operates at a temperature in a range from 4 degrees Kelvin to 77 degrees Kelvin.
 10. The system of claim 1 wherein each of the plurality of layers comprises a yttrium barium copper oxide (YBCO) material.
 11. The system of claim 1 wherein the layers are cooled using a coolant selected from the group comprising helium, neon and nitrogen.
 12. The system of claim 1 wherein the superconducting elements comprise a bulk superconducting material.
 13. The system of claim 1 wherein the magnet system has a homogeneous field distribution that varies less than 0.1 ppm for a 10 mm diameter spherical volume.
 14. The system of claim 1 further comprising a spacer positioned between a pair of the layers.
 15. The system of claim 14 wherein the spacer comprises a copper-silver alloy.
 16. The system of claim 14 wherein the spacer has a tensile strength of at least 500 MPa.
 17. The system of claim 1 further comprising a radio frequency (RF) device that emits a radio frequency signal within the cavity of the magnet system.
 18. The system of claim 1 wherein the superconducting layers have a tensile strength of at least 40 MPa.
 19. The system of claim 1 wherein the shim element further comprises a shim coil to adjust a magnetic field distribution within the cavity at room temperature.
 20. The system of claim 1 wherein the superconducting layers have a current density of at least 100 MA/m² at 77.3 degrees K.
 21. The system of claim 1 wherein the annular magnet has a field intensity of at least 5 T within a cavity.
 22. The system of claim 1 wherein a shim element is positioned between the superconducting annular layers.
 23. The system of claim 22 wherein the shim element comprises a coil in the cavity.
 24. The system of claim 22 wherein the shim element comprises a ferromagnetic tile.
 25. The system of claim 1 wherein the compound is selected from the group comprising DyBCO, NdBCO, HoBCO and YBCO.
 26. The system of claim 1 further comprising a processor to control operation of the system.
 27. The system of claim 26 wherein the processor is programmed to apply a stored pulse sequence and record magnetic resonance spectral data.
 28. The system of claim 1 wherein each of the plurality of layers comprises a plurality of thin film layers of semiconductor material having a thickness of less than 1 mm.
 29. A method of using a persistent mode magnet system for magnetic resonance spectroscopy comprising: energizing superconducting magnet layers to form a persistent mode magnetic field, the magnet layers comprising stacked annuli that surround a cavity for insertion of a material, the annuli including a rare-earth barium copper oxide compound having a tensile strength of at least 20 Mpa and a current density of at least 10 Ma/m² in a temperature range of 4k to 80k; subsequently inserting the magnet annuli into a magnetic resonance operating system generating a magnetic filed of at least 4.7 T; and obtaining a magnetic resonance spectrum of the material with the operating system.
 30. The method of claim 29 further comprising transporting the magnet from an energizing system to the operating system.
 31. The method of claim 30 further comprising cooling the magnet with a cooling system.
 32. The method of claim 31 further comprising cooling with a coolant selected from the group including liquid helium (e.g. T_(op)=4.2K), solid neon (e.g. T_(op)=<24 K), liquid neon (e.g. T_(op)=27 K), solid nitrogen (e.g. T_(op)=<40 K), liquid nitrogen (e.g. T_(op)=77 K).
 33. The method of claim 29 further comprising cooling the magnet during transport.
 34. The method of claim 33 further comprising transporting the annular magnet with a solid cryogen such as neon, nitrogen or liquid helium.
 35. The method of claim 29 further comprising providing an operating system including a conductive coil to adjust the magnetic field.
 36. The method of claim 29 further comprising using a shield during transport of the magnet.
 37. The method of claim 29 further comprising operating the superconducting layers having a current density of at least 100 MA/m² at 77.3 degrees Kelvin or lower.
 38. The method of claim 29 further comprising operating the magnetic resonance system at a temperature in a range of 4 K to 80 K.
 39. The method of claim 29 wherein the step of energizing the magnet further comprises: positioning an annular magnet in an energizing electromagnet; energizing the electromagnet to field cool the annular magnet; cooling the annular magnet with a cooling system to a selected temperature and discharging the electromagnet to trap a magnetic field; and inserting the annular magnet with the trapped field into the magnetic resonance system.
 40. The method of claim 29 further comprising programming a computer to control operation of the operating system.
 41. The method of claim 29 further comprising using an annular magnet including YBCO annuli.
 42. The method of claim 29 further comprising removing a bodily sample from a mammalian animal and using the magnetic resonance to obtain spectral data for the sample.
 43. The method of claim 29 further comprising inserting a pharmacological agent into the magnetic resonance system and obtaining a magnetic resonance spectrum of the agent.
 44. The method of claim 29 further comprising generating a magnetic field of at least 7 Tesla with the annular magnet.
 45. The method of claim 29 further comprising a magnetic field tuning device that comprises a shim coil.
 46. The method of claim 45 further comprising a spacer that has a tensile strength of at least 500 MPa. 